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Abstract 



D . 

Z/^ • The quark gluon plasma (QGP) at zero temperature and high baryon number is a system that 

may be present inside compact stars. It is quite possible that this cold QGP shares some relevant 
features with the hot QGP observed in heavy ion collisions, being also a strongly interacting system. 

jy^ ' In a previous work we have derived from the QCD Lagrangian an equation of state (EOS) for the 



cold QGP, which can be considered an improved version of the MIT bag- model EOS. Compared 



O ! to the latter, our EOS reaches higher values of the pressure at comparable baryon densities. This 

c/2 ' feature is due to perturbative corrections and also to nonperturbative effects. Here we apply 

this EOS to the study of neutron stars, discussing the absolute stability of quark matter and 

(N. ^ ^ 

^ . computing the mass-radius relation for self-bound (strange) stars. The maximum masses of the 

0^ ' sequences exceed two solar masses, in agreement with the recently measured values of the mass of 

o ■ 

^^ ' the pulsar PSR J1614-2230, and the corresponding radii of around 10-11 km. 
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I. INTRODUCTION 

In spite of the rapid progress in the field, the region of the QCD phase diagram with low 
temperature and high chemical potential is still not well understood. According to the cur- 
rent status, supported by different model calculations, there is a low temperature deconfined 
phase of quarks and gluons, the cold QGP, in which we may have color superconducting 
phases [l|. One of the open questions concerning the cold QGP is how free are quarks and 
gluons in this phase. In the simple picture based on the MIT bag-model quarks and gluons 
do not interact except when they hit the bag wall. In the opposite corner of the QCD phase 
diagram, i.e., at zero chemical potential and high temperature, the equivalent picture of a 
hot ideal gas of noninteracting quarks and gluons was dramatically changed after the series 
of experiments with relativistic heavy-ion collisions conducted at RHIC and now at LHC. 
In the new picture, quarks and gluons form a strongly interacting system in which nonper- 
turbative physics persist even after twice the critical temperature. In particular, the gluon 
condensates do not disappear very rapidly as previously expected. In the case of the cold 
QGP, experiments in laboratories cannot be performed directly, since compression of cold 
nuclear matter up to these tremendous densities cannot be achieved. However this compres- 
sion occurs presumably in the core of dense stars and the idea that we might find cold QGP 

hn 

in neutron stars has been around already for some decades [jt^- It is even conceivable that 
a whole star, not only its core, be made of quark matter 5| . This possibility will be explored 
in this work. 

The existence of a deconfined quark phase in the core of neutron stars [3|, Ij] depends 
crucially on the EOS. On the theoretical side there is still considerable freedom, since it is 
easy to calculate the mass and the radius of a star for a given EOS. Changing parameters in 
the proposed EOS one can arrive at rather different curves in the mass-radius diagram. On 
the observational side it is very difficult to obtain the mass and the radius of one single object. 
However, once this combined information is available it will provide a very strong constraint 
on the EOS of dense matter. The most recent data already impose some limits on the EOS 
parameters [6|-l8| . A previous analysis of the observational data from the neutron star EXO 
0748-676 presented in |9| concluded that most of the EOS are too soft and therefore unable 
to support the existence of neutron stars with a quark phase. In spite of this conclusion 



a 



being disputed [7|, new precise measurements of the pulsar PSR J1614-2230 carried out 



in [lO| yielded a mass of 1.97 ± O.OSMq for this object and led to the idea of interacting 



quarks if a core is present ll|] . Alternatively, a self-bound star, composed entirely of quark 

matter, could explain a massive pulsar if the pairing interactions and vacuum energy fall in 

11 
the right range [&\ , depending on the value of the radius which is still under discussion. It 

is then interesting to explore the existence of a self-bound deconfined quark phase made of 



noninteracting quarks [lOt|, as suggested in previous attempts [7|, 



3- 



In this paper we consider a quark star consisting of u, d and s quarks. Heavier quarks are 
not present in neutron stars [l3| . We shall further assume that the masses of the quarks are 
rriu = 5 MeV, rrid = 7 MeV, and rris = 150 MeV, complying with the generally accepted 
assumption of two light flavors and a heavier s quark. We first study the absolute stability 



parameter space of the EOS derived in IJ], which describes the quark gluon plasma at 



zero temperature. Stability requirements restrict the range of parameter values, which are 
subsequently used in the construction of the mass-radius diagram. 

This text is organized as follows. In Sec. II we briefly review the EOS for the cold QGP. 
In Sec. Ill we introduce the stability conditions and discuss its consequences. In Sec. IV we 
present the Tolman-Oppenheimer-Volkoff (TOV) equations for stellar structure calculations 
and their numerical solutions. Finally in Sec. V we present some comments and conclusions. 

II. EOS OF THE COLD QGP 

In [lj| EOS derivation started with the assumption that the gluon field can be decomposed 
into low ( "soft" ) and high ( "hard" ) momentum components. The expectation values of the 
soft fields were identified with the gluon condensates of dimension two and four, respectively. 
The former generates a dynamical mass, mc for the hard gluons, and the latter yields an 
analogue of the "bag constant" term in the energy density and pressure. Given the large 
number of quark sources, even in the weak coupling regime, the hard gluon fields are strong, 
the occupation numbers are large, and therefore these fields can be approximated by classical 
color fields. The effect of the condensates is to soften the EOS whereas the hard gluons 
significantly stiffens it, by increasing both the energy density and pressure. With these 
approximations it was possible to derive [14!] an analytical expression for the EOS, called 
here MFTQCD (Mean Field Theory of QCD). 

To proceed for the stellar conditions, we consider quarks u, d, s and electrons in chemical 



equilibrium maintained by the weak processes 15| : 



u + e -^ d + Ue 



u + e -^ S + Ue 



d ^ u + e + Ue 



S — )■ M + e + Ve 



and 



s + u ^ d + u. (1) 

As usual, the neutrinos are assumed to escape and do not contribute to the pressure and 
energy density. In chemical equilibrium we have 



A^d = yUs = /i and yUd + yUe = ^^■ 



The charge neutrality and baryon number conservation require 



gPu = gPd + -ps + Pe, 



and 



1/ 

PB = -^{Pu + Pd + Ps), 



(2) 



(3) 



(4) 



where pb is the total baryon density and pi is the density of quarks of flavor i {i = u, d, s) 
deflned by the corresponding Fermi momentum fcj given by 



pi 



27r2 



kf 



(note that we impose a local conservation of the charges). The electron density is 



(5) 



(6) 



where 7q and % are the quark and electron degeneracy factors given by 7q = 7e = 2 due 
to spin (the sum over color states was already performed). From ([T]) to (I?]) we flnd a set of 
four algebraic equations for Fermi momentum calculation for each particle: 



kj + k/ + ks^ = SttVb, 



^^M — kd I l^s ~r rCe , 



kd^ + rrid^ = ks^ + m^ 



(7) 



ku +mu^ + \Jke +me^ = yks +m 



for a flxed baryon density pB- The energy density is flnally given by 14 1 



/ 27g' 
VlGrriG^ 



Pb + Sqcd 



+ E 3 

i=u,d,s 



iQ \ h^\lh^ + "^i^ mi'^ki\Jki^ + m^^ nii 



27r2 



+ 



■In 



ki + \ki +mi 



Tfl 



+ 



7e I fee V/Ce + TWe^ m/key k^ + TUe'^ TTl 



27r2 [ 4 

and the pressure is 



-In 



nip 



ki + V fce' + me2 J + ^/n(me') [> , (8) 



P 



( 27^2 

1 16^;^; ^^ ~ "«^^ 



Pij' - Be, 



+ E 

i=u,d,s 



iQ jh yh +nii'^ 3miki\Jki + nii"^ 3m. 



27r2 



-In 



ki + yki +mi 



+ 



7g ( ke yke + nie"^ Smjk^yke + mg^ 3m, 



Gtt^ 



H — lnkf, + yke + m^^ 



3nii^ 
3ni ^ 



in{nii^) 



16 



ln{m 



e^) , 



(9) 
where me = 0.5 MeV" is the electron mass, ma is the dynamical gluon mass, and g is the 

coupling constant (a^ = 5'^/4vr) in QCD. Our analogue of the bag constant, called here 

Bqcd, is given by 

13qcd = ^K = {-/'""'F'P, (10) 

where (po is an energy scale associated with the energy density of the vacuum and with the 
gluon condensate 1J|. In ([8]) and iQ the summation over quark colors has already been 
performed. Throughout this work we employ the natural units h = 1, c = 1. 



III. STABILITY CONDITIONS FOR THE EOS 



We are interested in studying star models with stable strange quark matter. In this 
case, we have two stability conditions. The first one is that the energy per baryon of the 



deconfined phase 
matter defined in 



for P = and T = 0) is lower than the nonstrange infinite baryonic 



3, 



15| . Following these works we impose that: 



Ea = — < 934 MeV. 

PB 



(11) 



Since this condition must hold at the zero pressure point, from ([H]) and (Q we can numerically 
derive a relation between the bag constant Bqcd and the ratio ^ = g/mc- We solve iQ 
obtaining pB = pb{Bqcd, 0^ which is then inserted into IQ. The resulting expression is used 
to write the condition e^Bqcd, / Pb{Bqcd, = 934 MeV, which defines one " stability 
frontier ". This last equation is rewritten as ^ = C,{Bqcd), is plotted in Fig. 1 (solid line) 



and denoted by the 3-flavor line. Points in the Bqcd — C plane located on the right of 
the solid line are discarded since they do not satisfy flTTj) . The solid line, corresponding 
to the maximal value of E^ = 934 MeV, determines the maximum value of Bqcd — 75.7 
MeV/ fm?. The minimum value of Bqcd — 38 MeV/ fm? is determined by the second 
stability condition, which requires nonstrange quark matter in the bulk to have an energy 
per baryon higher than the one of nonstrange infinite baryonic matter. By imposing that 



Ea = — > 934 MeV 

PB 



(12) 



for a two flavor quark matter at ground state, we ensure that atomic nuclei do not dissolve 
into their constituent quarks. The constraint f lT2|) defines the dotted line in the Bqcd — i 
plane, denoted by the 2-flavor line in Fig. 1. Points located on the left of this line are 
excluded because they do not satisfy ( fT2l) . The region between the two lines in Fig. 1 
defines our stability window. 
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FIG. 1: Values of .^ = g/mc as a function of Bqcd for different values of the energy per baryon. 



The requirement of strange quark matter stability at finite pressure, in the interior of 
the stars, demands the introduction of another criterion. We shall assume that among the 
quark matter phase and the hadron phase, represented here by two hadronic models, the 
most stable is the one which has the highest pressure for the same value of the chemical 
potential. The curves p versus /x^ obtained with three EOS are shown in Fig. 2. As can 



be seen, our quark matter is more stable than the matter described by the hadronic models 
studied here. 
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ri'IG. 2: Pressure as a function of the chemical potential for the three EOS: MFTQCD ij] , Skyrme 



2l|, and Walecka 
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From Fig. 2 we can conclude that, at increasing chemical potential (and density), quark 
matter becomes more and more favored with respect to the hadronic matter studied here. 

We performed the causality check for Bqcv = 38 MeV/ fm~^ (close to the minimum 
value) and for (the maximum value) Bqcv = 75.7 MeV/fm~^. These two values of the bag 
constant define the stability range. Using these two values in Fig. 1 as entries to the dotted 
line and to the solid line, respectively, we can read in the vertical axis the two corresponding 
values of the variable ^, which are ^ = 0.007293 MeV~^ for Bqcv = 38 MeV/fm-^ and 
^ = 0.000657 MeV"^ for Bqcv = 75.7 MeV/ fm~^. Having fixed these parameters, we go 
back to (IH]) and (j9]) and, obtaining e and p for successive values of pb, we construct the EOS 
in the form p = p{e), plotted in Fig. 3. In this type of plot the slope is the speed of sound, 
which, due to causality, can not exceed the unity. This limit is shown by the full lines in the 
figure. 
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FIG. 3: EOS for the cold quark-gluon plasma. 

IV. NUMERICAL SOLUTIONS OF THE TOV EQUATION 

In order to describe the structure of a static, non-rotating compact star we solve the 
Einstein equations 16| : 

G^^ = -SttGT'^^ (13) 

for a spherical, isotropic, static, and general relativistic ideal fluids in hydrostatic equilib- 
rium. This particular solution of ( 1T3|) leads to the Tolman-Oppenheimer-Volkoff (TOV) 
equation for the pressure p(r): 



dp Ge{r)M{r) 

dr r^ 



1 + 



p(r) 



e(r) 



1 + 



47rr^j9(r) 



M{r-) 



1 



2GM{r) 



-\ -1 



(14) 



where G is the Newton gravitational constant. The enclosed mass M[r) of the compact star 
is given by the mass continuity equation: 

dM{r) 



dr 



Airr e(r). 



(15) 



Equations iHM and ( TT5l) express the balance between the gravitational force and the internal 
pressure acting on a shell of mass dM{r) and thickness dr. 

We solve numerically (fT^ and (TT5l) . which are coupled nonlinear equations for p{r) and 
M{r), to obtain the mass-radius diagram. The pressure and the energy density in (IT^ and 



8 



(fT5|) are given by the MFTQCD expressions ([9]) and ([8]), respectively. We take the central 
energy density to be e(r = 0) = ec and then we integrate out flT^ and flTSl) from r = up 
to r = R, where the pressure at the surface is zero: p{r = R) = 0. In Fig. 4 we show the 
mass-radius diagram for several values of Bqcv and ^ respecting the stability condition. 
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FIG. 4: Mass-radius diagram for several values of Bqcd and ^ allowed by the stability condition. 

Table |T] summarizes the values of mass, radius, and central energy density obtained for 
the several values of Bqcv shown in Fig. 4. At this point, the relationship of this parameter 
to other commonly employed quantities (i.e. the MIT bag constant) is difficult to assess, 
and the reasonable values obtained for the stability window to hold are quite encouraging. 



TABLE I: Bag, Maximum Mass, and Radius of the quark star. 



BQcv{MeV/fm 


3) ^(Mey-i) M{Mq) R{Km) 


62 


0.003658 2.56 11.97 


72 


0.002034 1.99 10.17 


75.7 


0.000657 1.82 9.69 



V. CONCLUSION 

In this paper we have apphed an EOS of the cold QGP to the study of compact stars. 
We note that when gluon interactions are switched off, we recover the standard MIT bag 
model EOS. The inclusion of gluon interactions generates more pressure and energy density, 
rendering the equation of state harder than the MIT bag model one and able to support 
stellar sequences with larger maximum masses. Indeed, our solutions of the TOV equations 
yield stars with two solar masses, in agreement with recent observations [lOj]. In the present 



paper we have improved a previous one 17|] in several aspects. The most important one was 
to introduce the requirement of stability, which strongly constrained the range of possible 
parameters. However, even after this strong restriction of parameter choice, we were still 
able to find stable quark stars with acceptable masses and radii. The latter is never too 
large {R < 12km), even for stellar sequences with maximum masses of ~ 2.5Mq, therefore 



determinations of radii underway 18|-|20| have the potential of constraining or even ruling 



out this type of theory in the near future. 
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